Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed January 9, 2013 (MN KTeX style file v2.2) 



Population synthesis of triple systems in the context of mergers of 
carbon-oxygen white dwarfs 

"on : A.S. Hamers 1 *, O.R. Pols 2 , J.S.W. Claeys 2 and G. Nelemans 2 ' 3 

1 Leiden Observatory, Leiden University, PO Box 9513, NL-2300 RA Leiden, The Netherlands 

2 Department of Astrophysics/IMAPP, Radboud University Nijmegen, PO Box 9010, NL-6500 GLNijmegen, The Netherlands 
- 'Institute for Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium 

Accepted for publication in MNRAS 2013 January 7. Received 2013 January 7; in original form 2012 November 21 

ABSTRACT 

Hierarchical triple systems are common among field stars yet their long-term evolution is 
poorly understood theoretically. In such systems Kozai cycles can be induced in the inner bi- 
nary system during which the inner orbit eccentricity and the inclination between both binary 
orbits vary periodically. These cycles, combined with tidal friction and gravitational wave 
emission, can significantly affect the inner binary evolution. To investigate these effects quan- 
titatively we perform a population synthesis study of triple systems and focus on evolutionary 
paths that lead to mergers of carbon-oxygen (CO) white dwarfs (WDs), which constitute an 
important candidate progenitor channel for type la supernovae (SNe la). We approach this 
problem by Monte Carlo sampling from observation-based distributions of systems within 
the primary mass range 1.0 - 6.5 M Q and inner orbit semi-major axes ci\ and eccentricities 
e\ satisfying a\ {\ - efj > 12 AU, i.e. non-interacting in the absence of a tertiary component. 
We evolve these systems by means of a newly developed algorithm that couples secular triple 
dynamics with an existing binary population synthesis code. We find that the tertiary signifi- 
cantly alters the inner binary evolution in about 24% of all sampled systems. In particular, we 
find several channels leading to CO WD mergers. Amongst these is a novel channel in which 
a collision occurs in wide inner binary systems as a result of extremely high eccentricities 
induced by Kozai cycles. With assumptions on which CO WD mergers lead to a SN la explo- 
sion we estimate the SNe la delay time distribution resulting from triples and compare to a 
binary population synthesis study and to observations. Although we find that our triple rate is 
low, we have determined a lower limit of the triple-induced SNe la rate and further study is 
needed that includes triples with initially tighter inner orbits. 
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1 INTRODUCTION 

Coeval stellar hierarchical triple systems, henceforth triple sys- 
tems, are kno wn to be quite common among stellar syst ems in 
the field (e.g. iTokovinin et al. I l2006l : iRaghavan et alll201fjh . Such 
systems consist of an inner binary pair orbited by a distant star, 
the tertiary. Due to secular gravitational effects high-amplitude 
eccentricity oscillations, known as Kozai or Kozai-Lidov cycles 
jKozdll962l :l Lido vl ll962h . can be induced in the inner binary sys- 
tem. These oscillations occur on timescales that are typically much 
longer than the orbital periods of both inner and outer binaries. 
One possible consequence of these eccentricity cycles is strong 
tidal friction and subsequent orbital shrinkage, a process which is 
known as Kozai cycles with tidal friction (KCTF). KCTF has been 
studied in the context of ( solar-mass) main sequence (MS) stars 
jMazeh & Shahaml 1 19791 ; lEggleton & Kiseleva-Eggletoi] 1200 ll : 
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20071) . In particular. Fabryckv & Tremaind "( 2007 ) have shown 
that KCTF is responsible for producing close MS binary systems 
with periods P] ~ 3 d, which is consistent with the observation 
that such systems are very likely (96%) orbited by a tertiary 
(Toko yinin et alj [2006). It remains to be seen, however, whether 
KCTF is still effective in higher-mass triples (individual masses 
m > 1.2 M ). This is because their constituents have radiative 
envelopes and are thus much less effective at dissipating tidal 
energy than their lower -mass counterparts, which have convective 
envelopes JZahn|[l977h . On the other hand, as such higher-mass 
stars greatly increase in size and develop convective envelopes 
during their red giant branch (RGB) and asymptotic giant branch 
(AGB) phases tidal friction is expected to become much more 
effective. Thus it is possible that KCTF is effective at significantly 
shrinking the inner binary system during the RGB/AGB phases, 
whereas this was not the case during the MS. 

Another possible consequence of high-amplitude eccentricity 
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cycles induced in (close) inner binary systems by a tertiary is the 
enhancement of the emission of gravitational waves (GWs), a pro- 
cess which i s strongly dependen t on the orbital eccentricity (e.g. 
|Peters|[l964 ). Thomp son] d201ll) has suggested that this mecha- 
nism can significantly reduce the GW merger time in close carbon- 
oxygen (CO) white dwarf (WD) binaries if they are orbited by a 
tertiary with the appropriate orbital parameters. This scenario im- 
plies an increase in the expected number of close CO WD systems 
that merge within a Hubble time, which is relevant for binary pop- 
ulation synthesis studies in which the CO WD merger channel is 
considered as an important candidate progeni t or channel for type 
la supernovae (SNe la) (e.g. iRuiter et alj20()A iToonen et al J2OI2I ; 
Claeys et al. in prep.). 

In this work we investigate both scenarios of KCTF during 
stellar evolution and the merging of CO WDs under the influence of 
a tertiary companion by means of a population synthesis of triples. 
We take into account long-term dynamical, stellar and binary evo- 
lution by means of a newly developed algorithm in which a code 
modeling the secular dynamical evolution of triples is coupled with 
the binary population synthesis code Binary.c (see Sect. 12.31 for 
details) to model the inner binary evolution. We restrict ourselves 
to intermediate mass triples with wide inner binary systems that 
would not interact in the absence of a tertiary component. For these 
triples we show in this paper that the presence of the tertiary com- 
ponent introduces to the inner binary system various novel evolu- 
tionary scenarios. In particular, we find an alternative pathway to 
SNe la as a result of collisions of WDs in wide inner binary sys- 
tems in which the eccentricity reaches extremely high values as a 
result of Kozai cycles. 

The organization of this paper is as follows. In Sect. [2] we de- 
scribe the triple evolution algorithm developed in this work and in 
Sect. [3] we provide details of two example systems that lead to CO 
WD mergers via two distinct triple channels. With this algorithm 
we perform a population synthesis study of which we discuss our 
methods in Sect.|4]and present the results in Sect. [5] For the various 
CO WD merger channels we find we estimate the expected triple- 
induced SNe la rates and compare these to the binary population 
synthesis study of Claeys et al. (in prep.) and to rates inferred from 
observations in Sect. [6] In Sec t.[7]we relate our results to the previ- 
ous work of lThompsorJ d20 1 lh and discuss the uncertainties of our 
treatment of triple evolution. We conclude in Sect. [8] 



2 HIERARCHICAL TRIPLE EVOLUTION ALGORITHM 

2.1 Secular hierarchical triple dynamics 

During Kozai cycles the mutual torque exerted by the inner and 
outer binary orbits causes an exchange of angular momentum be- 
tween the orbits while the or bital energies, and thus the s emi-major 
axes cij, are conserved (e.g. Mardling & Aarseth 2001). Here we 
use index j = 1, 2 to denote quantities pertaining to the inner and 
outer orbits respectively. As a result of the exchange of angular 
momentum the orbital eccentricities e ; and the total mutual incli- 
nation angle between both orbits i to( vary periodically on timescales 
P K » P2, where P , denotes orbital peri od. To lowest order in 01/02 
P K is given by dKinoshita & Nakaill999i) : 



P K = a 



P\ my + m 2 + m 3 
~P~i 



'"3 



(1- 



3/2 
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where m\ and m 2 denote the inner binary primary and secondary 
mass respectively and m 3 denotes the tertiary mass. The dimen- 
sionless quantity a is of order unity and depends weakly on the 



initial values of i tot , e\ and gi, where gj denotes the argument of 
periastron. Other quantities that vary periodically are gj (i.e. orbital 
precession, also known as apsidal motion) and the longitudes of the 
ascending nodes hj. 

In systems with large 01/02 the eccentric Kozai mech- 
anism may apply i n some cases dLithwick & Naoj 1201 lL 

Shappee & Thompson l2012h . This mechanism is associated with 
orbital "flips" (change of orbital inclination from retrograde to pro- 
grade and vice versa) with potentially very high inner orbit eccen- 
tricities. A quantity that measures the importance of this effect is 
the "octupole parameter" e oct defined by: 



m t - m 2 fl| 



( J 2 



mi + m 2 a 2 1 - ei 



(2) 



where typically logi (|e ocl |) 



-2 indicates that the eccen- 



tric Kozai mechanism cou ld be important dNaoz et al.l 1201 ll ; 
IShappee & Thompsonll2oT2T) . 

In physically realistic situations there exist sources of preces- 
sion in the inner orbit other than those due to secular three-body 
dynamics alone (in our study such effects in the outer orbit are neg- 
ligible because o 2 is always sufficiently large). We take into account 
the following additional contributions to precession: due to general 
relativistic effects (£i.gr)> due t0 tidal bulges induced in physically 
extended stars in relative close proximity (gi. t id e ) and due to intrin- 
sic spin rotation (gi, 10 , a , c ). The fact that g 1?GR , g l tide and gi. rotate are 
all positive implies that the associated sources of precession pro- 
mote precession in one particular direction and therefore suppress 
the resonance between gj and e\ during Kozai cycles. This phe- 
nomenon has been investigated in detail bv lBlaes et all d2002t) . 

Furthermore, whenever the distances between the components 
in triple systems become comparable to their radii tidal effects may 
become important. In our study such effects in the outer orbit can 
be neglected. However, in the inner binary system tidal effects may 
become important even for large semi-major axes with respect to 
the radii, oi //?,■ (where Rj denotes the inner binary stellar radius), if 
the inner orbital eccentricity is sufficiently high. In this case a tidal 
torqu e is induce d, either due to the equilibrium tide or dynamical 
tide dZahnll977h . that allows for the exchange between angular mo- 
menta of the spin and orbit in the inner binary system. This process 
continues until copl anarity, synchronization and circularization are 
achieved (Hut 1980). It is described quantitatively in terms of orbit- 
averaged equations for the orbital parameter time derivatives, o^tf, 
ei TF and fl/jr, where O, is the spin frequency of inner orbit star i 
and TF denotes tidal friction. 

Another effect we take into account is GW emission, which 
we model for the inner binary only. As a consequence of the emis- 
sion of GWs the inner binary is circularized and shrinks, which is 
described quantitatively by orbit-averaged equations for ej ow and 

Ol.GW- 



2.2 System of first-order differential equations 

We describe the secular Kozai cycles quantitatively in terms of 
coupled equations for g ; ,sTD and e^TD (where STD denotes sec- 
ular three-body dynamics), that follow after orbital averaging of 
the three-body Hamiltonian expressed in Delaunay's action-angle 
variables, applying Hamilton's equatio ns of motion and eliminating 
the longitudes of the as cending nodes lHarrington 1968; For d et al.l 
l200d : lNaoz etalj201lh . We employ equations for these four quan- 
tities accurate to octupole order, i.e. up to and including third order 
in oi/o 2 . In addition, an equation that expresses conservation of 
total angular momentum describes the change of ! to i. 
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We then combine all the physical effects discussed in Sect, l2.1l 
to form a system of eight first-order ordinary differential equations 
(ODEs) that reads: 



g2 
<2l 

<-'! 

fa 
8 

n 2 



= gl.STD + gl.GR + gl.tide + gl,rotate; 
= g2,STD; 

= fil.STD + ^l.TF + 6l,GW; 

- ^2,STD; 

= Ol.TF + «i,GW; 

= [G, (G, + G 2 8) + G 2 (G 2 + Gifl)] ; 

= "l.TFj 



(3) 



where S = cos(i' tol ) and G; denotes the orbital angular momen- 
tum. For expressio ns for the p r ecessi o n quantities #i,gr, guide and 
g i .ro tate we refer to | Blaes et al] (2002), iSmevers & Willems d200ll) 
and iFabrvckv & Trem aine (2007) respectively; forthe tidal evo- 
lution quantities £itf, <^i,tf and f2 lTF we refer to Hut and 
for the GW evolution quantities ei.ow and ai.ow we refer to Petersl 
( 1964). The main assumption in Eq.[3]is that all physical processes 
that are described are independent such that the individual time 
derivative terms can be added linearly. In addition, in the expres- 
sions for guide, £i,rotate, £i,n>> fli.TFi ^i,tf and 8 we assume copla- 
narity of spin and orbit at all times, even though Kozai cycles in 
principle affect the relative orientations between the spin and orbit 
angular momentum vectors and in turn a misal ignment of these vec - 
tors affects the Kozai cycles themselves (e.g. lCorreia et alj|201lh . 
We justify this assumption by noting that for the majority of sys- 
tems that we study the orbital angular momenta of both inner and 
outer orbits greatly exceed the spin angular momenta in magnitude, 
therefore the stellar spins cannot greatly affect the exchange of an- 
gular momentum between both orbits. 

The expression for 8 in Eq. [3] is derived by differentiating 
the relation that expresses conservation of total angular momen- 
tum, Gf Bt = Gj + G\ + 2G1G26, with respect to time and solv- 
ing for 8 in terms of Gj and Gj. Note that GW emission acts to 
change a\ and e\ and therefore affects the magnitude of Gi but 
not its direction. Although in principle tidal friction can also af- 
fect the direction of Gi if the stellar spins and orbit are not copla- 
nar, with our assumption of coplanarity this is not the case. This 

implies that, with Gj = Lj Jl - e^, where Lj = Lj(aj,mj) is the 
orbital angular momentum in case of a circular orbit and which 
is constant with regard to STD, the quantities Gj are given by 

G ; = -[ e ./(i- e 2)]. Gj ..^ sro . 



2.3 Coupling to binary algorithm 

The hierarchical triple evolution algorithm developed in this work 
couples the system of differential equations Eq. [3] to the existing 



rapid binary 
iHurlev et al.l 



J2004 120061 . 



wpul a tion s ynthesis algorithm Binary .c based on 
2000l [2002 ) with updates described in llzzard et al.l 



20091) and Claeys et al. (in prep.). This binary star 



algorithm includes stellar evolution and a variety of binary interac- 
tion processes such as mass loss due to stellar winds, mass trans- 
fer and common envelope (CE) evolution. We refer to Hurl ey et al.l 
j2002l) for further details of this binary population synthesis algo- 
rithm. The novel aspect of our approach is that the evolution of 
the inner binary system is taken into account through the coupling 
to Binary _c. We also take into account the evolution of the ter- 
tiary, which is treated as being isolated such that its mass may be 
computed from single stellar evolution. In this work we consider 
systems in which at all times the tertiary is sufficiently distant from 



the inner binary system such that it does not fill its Roche-lobe in 
its orbit around the inner binary system, nor it is involved in a CE 
phase with the inner binary system, a scenario whic h has been ex- 
plored in earlier studies (e.g. llben & Tutukovlll999r) . Therefore the 
latter processes are not modeled in the triple algorithm. 

In the ODE solver routine the system of differential equations 
Eq.|3]is solved for the duration of each time step of the Binary _c 
algorithm. In order to ensure proper convergence of the solutions 
of Eq. [3] which is quite stiff in nature, we use an ODE solver 
speci fically designed to integrate stiff ODEs JCohen & Hindmarshl 
1 19961) . We have checked the results of the ODE solver routine when 
not coupled with the Binary.c algorithm by computing the evolu- 
tion of the same hierarchi cal t riple systems as st udied in lFord et al.l 
< l2000l) . lBlaes et alJd2002l) and lNaoz et aflfeoill) and find very sim- 
ilar results for e,, ( t0 , and «i (where applicable) as function of time. 
We have also found near exact agreement of the time evolution of 
et and i tol when comparing to analytical solutions which exist in 
the quadrupole o rder approximation and w ith the further restriction 
that G 2 » Gi jKinoshita & Nakailfl999h . During each Binary.c 
time step quantities not included in the left hand side of Eq.[3]such 
as the inner binary masses and radii are assumed to be constant. 
Because the Binary _c algorithm uses adaptive time steps to match 
the rate of stellar and binary evolution it is always ensured that 
these quantities do not change significantly during each iteration of 
the triple evolution algorithm. In addition, if needed the triple algo- 
rithm decreases the Binary.c time step such that the inner orbital 
angular momentum Gi does not change by more than 2%. Thus, 
changes to the inner orbit semi-major axis and eccentricity due to 
secular three-body dynamics coupled with tidal friction and GW 
emission are communicated to the Binary.c algorithm at appro- 
priate times. 

The version of Binary.c used in this work enforces circular- 
ization at the onset of Roche-lobe overflow (RLOF) and CE evo- 
lution. A lthough this assum ption is generally justified for isolated 
binaries (Hurle y et al. I l2002l) . in triple systems high-amplitude ec- 
centricity cycles could be induced during the former phase, in par- 
ticular if the timescale of these cycles is comparable to or smaller 
than the timescale of mass transfer driven by RLOF. An accurate 
treatment of mass tr ansfer in eccentric or bits is beyond the scope 
of this work (see e.g. lSepinskv et al.l2007l) and we therefore choose 
to disable the ODE solver routine whenever mass transfer driven 
by RLOF or CE evolution ensues. Whenever the ODE solver rou- 
tine is active it is ensured that tidal evolution and GW emission in 
the inner binary system as calculated by Binary _c do not change 
the inner orbital parameters. This is necessary because in this case 
these processes are taken into account by the ODE solver routine 
through the equations for e t and «i in Eq.[3] In this manner the cou- 
pling of eccentricity cycles with tidal friction and GW emission in 
the inner orbit is calculated consistently. 

If the ratio /3 = a 2 ja\ is large the effects of STD are likely to 
be quenched by different processes in the inner binary system, in 
particular due to general relativistic precession in very tight inner 
binary systems, as discussed in Sect. |2.1| In such situations it is not 
necessary to solve Eq.[3]because the STD-related terms are negligi- 
ble and the inner binary evolution can therefore be fully handled by 
Binary.c. We therefore introduce a criterion that temporarily dis- 
ables the ODE solver routine and leaves the inner binary evolution 
to Binary _c in such cases. Specifically, the routine is disabled if 
P > 10y6 cl it.GR, where j8 cr i t ,GR is the value of f} such that the periods 
associated with precession due to STD (Eq.[T) and general relativity 
(approximately given by the reciprocal of gi.GR) are equal, which 
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we find is given by (omitting numerical factors of order unity): 

a/3 



2 \'/ 3 
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where c denotes the speed of light and G N denotes the gravitational 
constant. 

Outer orbital expansion due to mass loss in either inner or 
outer binary systems is taken into account with the assumption of 
fast, isotropic winds (i.e. Jeans mode) such t hat a 2 (m t + m 2 + m i) 
and e 2 remain constant dHuandfl956l . 1 1963b . With this relation a 
new value of a 2 is computed from the changes of m\, m 2 and m 3 
during each time step of Binary _c. 

Lastly, at each Binary _c time step the triple system is checked 
for dyn amical stability by means o f the stability criterion formu- 
lated by Mardli ng & Aarsethl d200 lh . including the ad hoc inclina- 
tion factor / = l-(0.3/Vr) i tot (with i tot expressed in radians). When- 
ever fj < /Scats where y8 clit is given by this stability criterion, the STD 
equations are no longer strictly valid we do not model the subse- 
quent evolution. The phase of dynamical instability would need to 
be computed with a three- body integrator code, which is beyond 
the scope of this work (see|Perets & K rattej2012l for a recent study 
of triple systems that become dynamically unstable). 



2.4 Fixed binary parameters and evaluation of various 
physical quantities 

We detail some of the fixed binary evolution parameters used in this 
work. We choose quasi-solar metallicity, i.e. Z = 0.02. The initial 
spin periods of t he inner orbit comp onents are computed from a 
formula given by Hur ley et al which is fitted from data of 

the rotational speed of MS stars of Land d 19921) . The CE parame- 
ter q-ce describes the efficiency of the conversion of orbital energy 
into binding energy with which to shed the envelope and it is set 
to the canonical value of a C E = 1. In this CE description an addi- 
tional parameter Aqe is required, which is a dimensionless measure 
of the relative density distribution within the envelope; here it is 
determined by means of functional fits as detailed in Claeys et al. 
(in prep.). All other parameters intrinsic to the Binary _c algorithm 
are also set to identical values as in Claeys et al. (in prep.). 

In addition, various physical quantities are required to evalu- 
ate the right-hand sides in Eq.[3]and in the remainder of this sec- 
tion we describe how these quantities are obtained in our triple al- 
gorithm. The tidal evolution equations (i.e. expressions for aijp, 
ei,TF and £1i,tf) contain a parameter that captures the e fficiency o f 
tidal friction and depends on the structure of the star (Hutl ll98lh . 
This parameter is the ratio k ami /Tj, where k. imi is the classical ap- 
sidal motion constant of inner orbit star i and T, a tidal friction 
timescale. We compute fca m.,/7*, according to the same prescription 
that is used in Binary.c (Hurley et al. 20o3). In this prescription 
a distinction is made between convective, radiative and degenerate 
envelope s and k. dmJ /T ,- for t h ese three cas es i s computed based on 
results of R asio et al ] dl996)) . IZahnl d 19771) and lCampbelH d 19841) re- 
spectively. Furthermore, the gyration radii r 2 . are required for the 
evaluation of £%_tf and are computed with a prescription in which 
the stars are split into two parts, cons istin g of the core and the enve- 
lope, detailed in lHurlev et ail d2000h and lHurlev eTaLl d2002h . The 



quantities A; am ,/r, and rz t are therefore calculated from precisely 
the same prescription as in Binary.c, which is required for a con- 
sistent treatment. 

The classical apsidal motion constant £ am> , (i.e. not the ra- 
tio k am j/Tj), which is requ i red for evaluation of guide a nd gi, rotate 
dSmevers & Willemsll200ll : iFabrvckv & Tremainej|2007l) . is com- 
puted as a function of mass and time relative to the ZAM S from 
detaile d stellar models for metallicity Z = 0.02 calculated bv lClaretl 
(2004). The run of k. im , with time is determined by means of linear 
relations scaled to the time spent during each evolutionary stage 
for all tabulated masses in Claret (2004). Subsequently, A: amj is 
calculated for arbitrary mass by means of linear interpolation be- 
tween adjacent mass values. Furthermore, for low-mass stars, i.e. 
nij < 0.7 M , the classical apsidal motion constant is taken to be 
k am .i = 0.14 corresponding to n = 3/2 polytropes as an estimate 
for these (nearly) fully convective stars. For helium (He) MS stars 
and WDs k am j is calculated as fu nction of m ass by means of in- 
terpolations of tabulated data from lVilaldl977l) . Lastly, for neutron 
stars k am j is calc ulated as function of mass and radius from the ex- 
pression given in Hinderer (2008) and for black holes k dm j = as 
appropriate for non-spinning black holes. 



3 TRIPLE EVOLUTION EXAMPLES 

To demonstrate the triple algorithm presented in Sect.|2]we describe 
in detail the evolution of two triple systems in which an inner bi- 
nary CO WD merger occurs in either a tight circular orbit or a wide 
and highly eccentric orbit. The initial conditions for these two sys- 
tems are selected from our first triple population synthesis sample 
(TSMl.cf. Sect.|4~2l 



3.1 Merger in a tight circular orbit 



We consider a triple system with initial parameters mi = 3.95 M , 
m 2 = 3.03 M Q , m 3 = 2.73 M Q , a y = 19.7 AU, a 2 = 636.1 AU, 
ei = 0.23, e 2 = 0.82, i m = 116.0°, gi = 28.5° and g 2 = 249.6°. 
Fig.rjj shows fli, ei, P = 02/01 , the tidal strength quantity k^j/Tj 
dHutf[l98ll ; l"Hurlev et alj|2002h and i tol as a function of time. The 
evolution is shown with both the triple algorithm enabled (triple 
case; solid line) and disabled (binary case; dashed line), leaving 
all inner binary evolution to Binary.c. The latter binary case is 
equivalent to the situation in which no tertiary is present. 

In the binary case the semi-major axis increases at the two mo- 
ments when the binary stars evolve from AGB stars to CO WDs (at 
t x 0.2 Gyr and / x 0.5 Gyr respectively) as a consequence of the 
associated wind mass loss. Here (i.e. in Binary.c) it is assumed 
that the wind mass loss is fast and isotropic such that £?i(mi + m 2 ) 
remains approximately constant. When either star in the inner bi- 
nary system enters the RGB or AGB phase it develops a convec- 
tive envelope and its radius increases significantly. The change of 
the envelope structure is reflected by a substantial increase of the 
tidal strength quantities k-^j/Tj as is shown in Fig. Q] (primary 
AGB phase). However, despite the large k am i /Ti and primary ra- 
dius there is no significant tidal friction during the primary AGB 
phase because of the large a\ and small (and constant) e\. Conse- 
quently, in the binary case the CO WDs end their evolution in a 
wide and non-interacting binary with a\ « 84 AU. 

In the triple case high-amplitude Kozai cycles are induced dur- 
ing the MS in the inner orbit by the tertiary with maxima of ei, max ~ 
0.9. Octupole order terms are important because |e oct | a 10" 2 . Nev- 
ertheless, during the MS the eccentricity is not high enough to drive 
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Figure 1. Several quantities of interest in the evolution of the first example system (Sect. 13. U . Shown as a function of time are a i , ei , /? = ailai , k- dm j/Tj 
and ( to t. In the plots in the first column of the first two rows the entire evolution is shown; the other plots correspond to the end of the primary core helium 
burning phase, the primary AGB phase and start of the primary CO WD phase. Solid line: triple case; dashed line: binary case (where applicable). In the plot 
of k^mj/Tj only the triple case is shown (the binary case is very similar); in this plot the solid line applies to the primary (i = 1) and the dashed line to the 
secondary (;' = 2). Note that the evolution is not fully sampled in all plots, causing several kinks (in particular in the plots for e\ and i tot ) - more detailed 
calculations are performed internally in the ODE solver routine but are not shown here (cf. Sect. |2.3t . 



significant tidal friction because the tidal strength quantity is very 
small (k dm j/Tj ~ 10~ 18 s~') as a consequence of the radiative en- 
velopes in the MS stars. During the primary RGB phase the pri- 
mary possesses a convective envelope (k :imi /T 1 ~ 10~ s s~') but 
e\ is not high enough to trigger significant tidal friction. However, 
tidal friction does become effective during the primary AGB phase 
starting at f « 220.6 Myr and circularizes the inner orbit during 
the time span of five Kozai cycles, where significant orbital shrink- 
age occurs at eccentricity maxima. Note that during this phase the 



increase in the effectiveness of orbital shrinkage is due to the ex- 
pansion of the primary from a radius of R[ « 49 R Q to Ri a 497 R Q 
between f = 220 Myr and / = 222.5 Myr. Consequently «[ is re- 
duced to a i ss 7 AU and the orbit is completely circularized. Note 
that for complete circularization to occur the duration of the phase 
in which k. dm \IT\ is substantial must be sufficiently long compared 
to the Kozai period P K (Eq.QJ, which is the case for this example 
system. In other cases of triple systems, however, P K can be much 
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Figure 2. Quantities of interest for the second example system in which an inner binary CO WD merger occurs in a wide and highly eccentric orbit. 



longer than the duration of the RGB/AGB phases, thus avoiding 
strong tidal friction even if the eccentricity maxima are high. 

Shortly after KCTF during the primary AGB phase the pri- 
mary has swelled up to R\ » 61 1 /?o, fills its Roche-lobe and 
invokes CE evolution, shrinking the orbit further to «i ss 1 AU 
(r « 222.6 Myr). During this process the ratio /? increases substan- 
tially to /3 ~ 1321, which is large enough for general relativistic 
precession to fully quench any subsequent Kozai cycles. Conse- 
quently, the inclination angle ( tot is frozen to i tot x 128°. The pri- 
mary emerges from the CE as a CO WD and the secondary as a 
MS star. At t x 0.5 Gyr the secondary evolves to an AGB star it- 
self, invoking a second CE phase. A close CO WD binary with 
fl[ w 0.01 AU emerges from the CE. Due to GW emission the 
system subsequently merges at t « 13.2 Gyr in a circular orbit. 
The combined CO WD mass, 1.45 M Q , exceeds the Chandrasekhar 
mass M C h ~ 1.4 M , therefore the merger could potentially lead to 
a SN la explosion (cf. Sect. 16.2b . 



3.2 Merger in a wide eccentric orbit 

The second example system has initial parameters m\ = 3.26 M G , 
m 2 = 2.00 M Q , m 3 = 1.39 M ,ai = 716.6 AU, fl 2 = 2.012- 10 4 AU, 
ei = 0.78, e 2 = 0.64, i tot = 93.4°, gi = 150.0° and g 2 = 151.2°. The 
evolution of a t , e\,fi = a 2 la\ and i tot for this system as function of 
time is shown in Fig. ff] The main difference in initial parameters 
in this example system compared to the previous example is that a\ 
and a 2 are larger. Their ratio /? is comparable, however. Because of 
the larger a 2 the Kozai period P K « 100 Myr is much longer than 
in the first example (P K w 0.3 Myr). Consequently, it is much more 



unlikely that the primary RGB and AGB phases (which occur at 
t x 0.3 Gyr and t x 0.4 Gyr respectively) coincide with eccentricity 
maxima. Such a coincidence is indeed not the case in the second 
example system. 

As the primary and secondary evolve to CO WDs without in- 
teracting the associated wind mass loss widens the orbit to a\ x 
2600 AU. Although this inner binary wind mass loss increases a 2 
as well, the ratio j} = a 2 ja x decreases with the assumption of fast 
isotropic winds such that a\(m\ +m 2 ) and a 2 (m\ +m 2 + mi) are con- 
stant (in this example wj also remains constant). Consequently, the 
ratio p decreases to/? a 18 after the secondary has evolved to a CO 
WD. As the tertiary evolves to a CO WD and subsequently loses 
mass at t as 1.5 Gyr a 2 increases while a\ stays constant, thus in- 
creasing p again. The latter process also increases the Kozai period. 
After I » 1.5 Gyr the maximum values of cos(i tot ), which are still 
negative (retrograde orbit), slowly increase and the corresponding 
eccentricity maxima increase as well. As cos(i tol ) approaches the 
eccentricity maxima reach very high values until at t x 7.0 Gyr 
the orbits become prograde, which is associated with e l reaching a 
value of 1 - ei ss 10~ 74 . This is high enough to trigger an orbital 
collision between the two white dwarfs, even though a\ w 2600 AU 
is very large. Such a collision could potentially result in a SN la ex- 
plosion as in the previous example (cf. Sect. 16.2) , but note that the 
nature of the merger prior to a possible SNe la event is very differ- 
ent. 

One might wonder whether in this system tidal friction is sig- 
nificant during close periastron passages before the collision oc- 
curs. For this reason we have recalculated this system with the tidal 
strength quantities k.^j/Tj for radiative and degenerate damping 
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multiplied by ad hoc factors of 10 3 and 10 6 to artificially increase 
the effectiveness of tidal friction. Even with these multiplications 
we find that the outcome is identical, showing that tidal friction is 
not important in this system. We have performed similar tests for all 
systems in the population synthesis study and have found no effect 
(see Sect.lTIt. 



4 POPULATION SYNTHESIS: METHODS 

Before presenting the results of our population synthesis study 
(Sect.O we discuss here our selection criteria and sampling meth- 
ods which are essential ingredients for such a study. 



4.1 Selection criteria 

Because we are mainly interested in inner binary mergers involving 
a CO WD we take an upper limit of the inner binary primary mass 
rti\ of m\ „ = 6.5 M G , which approximately captures the boundary 
between forming a CO WD and an ONe WD in the case of single 
star evolution. The lower limit of mi is set to m u = 1.0 M Q for 
any lower mass will not produce a significant number of CO WDs 
within a Hubble time for a metallicity of Z = 0.02. In addition we 
focus on triple systems with wide inner binary systems such that 
the inner binary stars do not interact during their evolution in the 
absence of a tertiary. Taking into account possible tidal friction this 
can be achieved by requiring that the initial inner orbit semi-latus 
rectum /[ = ai (l — efj satisfies li > = 12 AU. 

One possible method to obtain the initial triple parameters is 
to directly sample data from observed triple systems. Fig. [3] (top) 
shows the observed sample of 725 triple systems of iTokovininl 
d2008h in the (Pi , /^-diagram. A clear paucity of systems with in- 
ner orbital periods between log 10 (P[/d) = 2 and log^Pj/d) = 4 is 
present, which is likely because of the difficulty of detecting such 
systems in the regime between spectroscopic and visual binaries. 
Our selection criteria of inner binary systems with 1.0 < mJM G < 
6.5 and /[ > 12 AU decrease the number of systems from 725 to 
165, where a value of log 10 (f i/d) = 4 is taken to represent the 
dividing value l t J = 12 AU. The number of remaining systems 
(indicated in the top figure of Fig. [3] by dots) is not large enough 
for direct sampling, therefore we choose to utilize the method of 
Monte Carlo sampling from observation-based distributions. In or- 
der to gain information on the uncertainties of our results we use 
two distinct sampling methods that are described below. 



4.2 Sampling methods 

4.2.1 TSM1 

The first triple sampling method (TSM1) is based on the supposi- 
tion that a hierarchical triple system is composed of two uncorre- 
lated binary systems. The main advantage of this approach is that 
the statistics of binary parameters are known with greater certainty 
than those of trip le systems. A primar y mass 1.0 < mJM e < 6.5 
is sampled from a [Kroupa et alj (1993) initial mass function (IMF), 
i.e. dN/dmi cc m~ 2 - 70 . Subsequently, two mass ratios, q\ = m 2 /mi 
and q 2 = m i /(m l + m 2 ), are sampled independently from a uni- 
form distribution such that < qj < I, m 2 > 0.01 M Q and 
m-i > 0.01 M e , consistent with Claeys et al. (in prep.). The lower 
limit on the stellar masses is in accordance with the limit found by 
iKouwenhoven et al] (2007). The secondary and tertiary masses are 
subsequently computed from these sampled mass ratios. We sample 
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Figure 3. Inner (Pi) and o uter (P 2 ) orbital p eriods of triple systems of 
the observed triple sample of Tokovinin (2008) (top) and the two synthetic 
triple samples described in Sect. I4.2l (middle and bottom). In the observa- 
tional sample systems satisfying our selection criteria are denoted with large 
(red) dots (•). In the synthetic samples the large dots (•) denote the systems 
that undergo an early MS inner binary merger or destabilization and are 
therefore not likely to be observed. The line corresponds to P 2 /Pi = 5 and 
represents the typical critical value for dynamical stability. 



ei and e 2 from a thermal distribution diV/de, = 2e j as is appropriate 
for Pj > 10 3 d binaries dKroupa & BurkerfeOOlh . Two semi-major 
axes are then sampled from a distribution that is flat in log 10 (ay), the 
smaller of which is designated a t and the larger of which is desig- 
nated a 2 . The lower and upper limits a t and a u in the distribution of 
the semi-major a xes are a; = 5R and a u = 5 • 10 6 P Q for b° m inner 
and outer orbits dKouwenhoven et al]|2007l) . From these systems 
we reject those that do not satisfy l\ = a-i (1 - e 2 ) > l\j = 12 AU. 
Lastly, we sample the initial orbital inclination angle i tol from a 
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Tokovinin TSM1 TSM2 

0.25 

Mean SD Mean SD Mean SD 



mi/M e 


2.36 


1.23 


1.85 


1.01 


2.17 


1.24 


m 2 /M 


1.69 


0.94 


0.93 


0.79 


1.09 


0.93 




1.53 


1.39 


1.36 


1.23 


1.59 


1.43 


log 10 (Pi/d) 


5.13 


0.73 


5.47 


0.81 


5.77 


1.29 


log 10 (P 2 /d) 


7.04 


0.76 


7.87 


0.75 


8.28 


1.61 



Table 1. Mean and standard deviation (SD) of masses and orbital periods 
of the lTokovininl |2008) sample satisfying the selection criteria (Sect. |4~T1 
and of both sampled triple populations TSM1 and TSM2. 

distribution that is uniform in cos(z tot ) with -1 < cos(; tot ) < 1 
and the initial arguments of periastron of both orbits, g\ and g2, 
from a uniform distribution with < gj < 2n. From the triple sys- 
tems obtained in this manner we subsequently reject the systems 
that are not dynamically s table based on the stability criterion of 
Mard ling & Aarsethl (200 1 ), consistent with our triple evolution al- 
gorithm (cf. Sect.|2"3t. 

4.2.2 TSM2 

In the second triple sampling metho d (TSM2) we use the mul- 
tiple system recipe developed by lEggletonl (2009). This recipe 
is designed to reproduce the properties of a set of 4558 stellar 
systems with Hipparc o s m agnitude brighter than 6 collected by 
lEggleton & Tokovininl J2008I) . It is followed until the second bi- 
furcation, effectively restricting to multiple systems that are triple. 
The parameters given by the recipe are m\, mi, m^, Pi and P 2 . The 
primary mass distribution in this recipe is designed to resemble the 
Salpeter IMF (dN/dm l oc m^ 235 ) at large masses and turns over 
for masses below 0.3 M Q . The mass ratio distribution is approxi- 
mately flat and the P2 distri bution has a broad peak around 10 5 d. 
We refer to Eggleton (2009) for further details. The corresponding 
semi-major axes are computed from Kepler's third law. The pa- 
rameters not prescribed by the recipe are i tot , ej and gj, which are 
sampled from the same distributions as in TSM1 above. Similarly 
to the method in TSM1 systems are rejected if l\ < lu = 12 AU and 
if they do not satisfy the stability criterion of Mardlin g & Aarsethl 
( 120011) . 

4.2.3 Properties and comparisons 

Here we elaborate on some of the properties of the sampled and ob- 
served populations. Fig.[3](middle and bottom) shows the (Pi, P2)- 
diagram for a small sample (360 systems) of the TSM1 and TSM2 
populations. The upper boundaries of Pi and P2 in both populations 
are quite different. In TSM1 Pi and P2 are limited by a\ and a„, 
whereas in TSM2 such sharp boundaries do not exist. In TSM2 the 
upper boundary of P2 is dependent on Pi because in this sampling 
method the largest ratio P 2 /Pi is given by P 2 /Pi = 2 ■ 10 6 . These 
wide outer binary systems are very rare, however. Fig. |3](middle 
and bottom) also shows the systems that according to our triple al- 
gorithm experience a merger in the inner binary system or a de sta- 
bilization early on the MS (cf. Sects. |5TTTl~l and [5~3"l l. Such systems 
are not likely part of an observed sample of triples because of their 
short lifetimes. Because of their small number, removing these sys- 
tems does not significantly affect the distributions of masses and 
orbital periods (the mean values are affected by typically 1%). 

In Table[T]the distributions of the masses and orbital periods of 




io gl0 (iypj) 

Figure 4. Probability density functions (PDFs) of P2/P1 for the Tokovinin 
sample with our selection criteria applied (solid line), the TSM1 sample 
(dashed line) and the TSM2 sample (dotted line). For the Tokovinin sample, 
errors based on Poisson statistics are indicated. 

the Tokovinin sample with our selection criteria applied are com- 
pared to the two sampled populations. Because of the steep slope 
of -2.70 in the TSM1 IMF the mean masses in TSM1 are small 
compared to those in TSM2, for which a Salpeter-like slope (-2.35) 
applies for mi > 1.0 M Q . Compared to the Tokovinin sample the 
masses of TSM1 and TSM2 appear to be on the low side, but it is 
important to note that the Tokovinin sample is magnitude-limited 
rather than volume-limited, therefore it is strongly biased towards 
more massive systems. The distributions of the orbital periods of 
TSM1 and TSM2 are similar, the orbital periods being on average 
somewhat higher in TSM2 compared to those in TSM1. The orbital 
periods of the Tokovinin sample, notably the outer orbital periods, 
are smaller than those in the sampled populations. This is also likely 
affected by observational bias, because very long orbital periods on 
the order of 10 7 days or longer are generally very difficult to mea- 
sure. 

Lastly, Fig. [4] compares the distributions of the ratio P2/P1 
for the three populations. For TSM1 and TSM2 these distributions 
are similar. The distribution of P2/P1 of the Tokovinin sample is 
peaked towards smaller ratios, although here observational bias for 
large period ratios may be important. Note also the large error bars 
in the Tokovinin sample as a consequence of the rather small num- 
ber of systems in the Tokovinin sample that satisfy our selection 
criteria (165). 

4.2.4 Sampling procedure 

We sample = 2 ■ 10 6 triple systems for both populations TSM1 
and TSM2. Using the triple algorithm described in Sect. [2] these 
systems are evolved from t = 0, with all three components starting 
as ZAMS stars, to a Hubble time t H = 13.7 Gyr. In order to retain 
sufficient resolution in higher-mass systems (2.0 < mi/M Q < 6.5) 
we split both TSM1 and TSM2 into two parts of 1 • 10 6 systems, one 
with 1.0 < mi/M Q < 2.0 and one with 2.0 < mi/M Q < 6.5. In the 
results presented in Sect[5]systems in each part are given appropri- 
ate weights determined by the mass function of mi to account for 
the fact that the actual number of systems occurring in the 1.0 < 
mi/M < 2.0 range is larger than that in the 2.0 < rai/M < 6.5 
range (see AppendixfA]for details). Several key events in the evolu- 
tion are kept track of, most notably the onset of significant KCTF, 
CE evolution, a merger in the inner binary system and a desta- 
bilization of the triple system. We thus obtain a catalogue of the 
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Description 


o 


MS (m < 7 


1 


MS Im > 7 Af^l 


2 


Hcrtzsprung (HG) 


i. 

j 


red giiint branch (RGB) 


4 


r*rvrf nflmm niirnino li TTpRi 


5 


early asymptotic giant branch (EAGB) 


6 


thermally pulsing AGB (TPAGB) 


7 


naked helium star MS (He MS) 


8 


naked helium star Hertzsprang gap (He HG) 


9 


naked helium star giant branch (He GB) 


10 


helium white dwarf (He WD) 


11 


carbon-oxygen white dwarf (CO WD) 


12 


oxygen-neon white dwarf (ONe WD) 


13 


neutron star (NS) 


14 


black hole (BH) 


15 


massless remnant 



Table 2. Descriptio n of the different stellar types used in Fig. [5] Identical to 
the types used in lHurlev et~ai] J2002h . 

evolutionary outcomes where we distinguish between three main 
channels: inner binary mergers, no inner binary mergers and triple 
destabilizations. In a very small number of systems (55 and 44 for 
TSM1 and TSM2 respectively) the calculation of the evolution can- 
not be completed due to convergence errors in the ODE solver rou- 
tine (cf. Sect.0; these systems are excluded from the results below. 
In these cases the errors occur just prior to a MS merger or RGB + 
MS merger. Because the affected systems constitute a tiny fraction 
of all systems for which a MS merger or RGB + MS merger applies 
this does not affect our conclusions. 



5 POPULATION SYNTHESIS: RESULTS 

We present the main channels that we find in our population synthe- 
sis study in Table[3] where the first two columns show the probabil- 
ities P (in per cent) of the main channels and several sub channels 
of both sampled populations TSM1 and TSM2. In addition to this 
table, Fig. [5] shows more detailed information on the likelihood of 
the different combinations of inner binary mergers where the stellar 
types used are defined in Table|2] A post-MS merger is defined as a 
merger between a post-MS (non-compact object) primary and any 
secondary, whereas a compact object merger is defined as a merger 
between a compact object primary and any secondary. In addition 
to the probabilities of the channels information on the occurrence 
of KCTF, CE evolution and the nature of the merger (i.e. occur- 
ring in a circular or highly eccentric orbit) is given by means of 
the quantities /kctf, /ce and /ec- We define /kctf as the fraction 
of systems in each channel in which I^lstdI > 10~ 18 s~' (such that 
there is a significant change in e\ due to STD in a Hubble time) and 
0.1 |ei,TFl < I^i.stdI < I^i.tfI (i-e. ei,sTD and eijp are of comparable 
order of magnitude). This can happen at various points in the evolu- 
tion and is strongly influenced by the envelope structure of the inner 
binary stars because this determines the effectiveness with which 
tidal energy can be dissipated. In particular if any inner binary star 
possesses a convective envelope (most notably low-mass MS, HG, 
RGB and AGB stars) strong tidal friction is likely if e t is also high 
due to Kozai cycles. In Table [3] a distinction is made between the 
types of the star in which tidal energy is dissipated for the duration 
of KCTF. Furthermore, f C E is defined as the fraction of systems in 
each channel in which CE evolution is at some point invoked. This 



includes both CE evolution triggered in highly eccentric orbits and 
CE evolution triggered by RLOF (in circular orbits). The fraction 
fcE includes the possibility for multiple phases of CE evolution. 
Lastly, the quantity /ec is defined as the fraction of systems in each 
channel in which a merger occurs due to an eccentric collision. We 
define such a collision to occur if the inner binary periastron dis- 
tance is smaller than or equal to the sum of the inner binary radii, 
i.e. iffli(l-ei)<i?i+i?2 and when in the corresponding circular 
case there would not be a collision, i.e. a { > R { + R 2 . 

An important result of our simulations shown in Table|3]is that 
in as 8% of all systems an inner binary merger occurs, in ss 6% of 
all systems no inner binary merger occurs but the orbit is shrunk 
significantly and that ss 10% of all systems become dynamically 
unstable at some point in the evolution. This means that in « 24% 
of all systems the tertiary significantly alters the inner binaiy evo- 
lution whereas in the absence of the tertiary there would not be 
such interaction. Taking into account that systems in which the in- 
ner binary components merge during the MS (~ 4%) or in which a 
dynamical instability occurs during the MS (~ 4%) are not likely 
to be observed (cf. Sect. I5.lt this estimate is reduced by » 8%. 
This implies that the tertiary significantly alters the inner binary 
evolution in as 16%/(1 - 0.08) ~ 17% of the "observable" triples, 
i.e. triples in which the hierarchical structure is not disrupted in an 
early stage in the evolution. 

A second important feature evident from Table [3] is that the 
results for TSM1 and TSM2 are generally similar, even though they 
have been sampled by different methods. To provide insight into 
the dependence of the probabilities of the main channels on the 
initial parameters we show in Fig. |6]these probabilities (expressed 
as fractions) as function of the initial parameters q\ = m2jm\, a\, 
P and ( tot . We will refer to this figure and Table|3]when discussing 
the three main channels in more detail below. 



5.1 Inner binary mergers 

As a consequence of high-amplitude eccentricity cycles induced 
in the inner orbit the components in the inner binary may at some 
point merge, either directly because of an orbital collision (included 
in Table[3]in the fraction / EC ) or indirectly as a result of strong tidal 
friction induced by high eccentricity, possibly invoking CE evolu- 
tion. In the latter case the inner binary can either merge during the 
CE or survive the CE in a tight orbit and (possibly after another CE 
invoked by the secondary) subsequently merge as a result of orbital 
energy loss due to GW emission. We separate further discussion 
into mergers occurring on the MS, after the MS and mergers in- 
volving a primary compact object. 

5.1.1 MS mergers 

As Table [3] demonstrates most mergers occur on the MS and for 
almost all of the latter the merger occurs through an eccentric col- 
lision. In the latter case -4 < log 10 (l - e\) < -5 just prior to 
merger. In the Binary.c algorithm it is assumed that the result of 
the merger is a single MS star with mass m\ + mi\ depending on 
the latter value the final remnant is either a CO WD (most cases), 
an ONe WD or a NS (following a core-collapse supernova). We 
find that octupole terms in the STD are very important for MS 
mergers, which is reflected by the large values of |e oct | (cf. Eq. O 
prior to merger for these systems; the latter quantity is narrowly 
distributed around |e ocl | « 10 _1 5 . This can be attributed to a com- 
bination of initially small semi-major axis ratios p = a 2 la\, small 
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Figure 5. Left: the number of systems in TSM1 that experience an inner binary merger for each occurring combination of stellar types just prior to merging, 
normalized to the total number of merger systems (that constitute as 8% of all sampled triple systems). Right: the number of systems in TSM1 in which 
the merger occurs d ue to an eccentric collision, again normalized to the total number of merger systems. The stellar types are identical to the ones used in 
iHurlev et all 120021) and are repeated for convenience in Table|2] 



inner binary mass ratios qi = m^jmx and high outer orbit eccen- 
tricities e?- Fig. [6] indeed shows that the fraction of MS mergers, 
/ MSmergc , decreases significantly with increasing q t and is peaked 
towards small p around 0.5 < log 10 (/?) < 2.0. The MS mergers 
are also characterized by high initial inclination angles (cf. Fig. 
[SJ. Note that a maximum in /ms merge occurs at inclinations slightly 
above i lot = 90°, i.e. fo r retro grade orbits, which is consistent with 
Shapp ee & Thom pson (20l2) who find an asymmetry in cos(7 tot ) 
towards negative values for the probability of the eccentric Kozai 
mechanism to occur. With regards to the dependence of /ms merge 
on flj it might be expected that MS mergers likely occur for small 
initial a\ as a tighter orbit makes an eccentric collision more likely. 
However, we find that MS mergers on average have large initial 
with a maximum of /ms merge around a\ « 10 28 AU ~ 6 • 10 2 AU 
(cf. Fig- [6j. This implies that the effect of increasing the eccentric- 
ity maxima with decreasing /? (and thus increasing a{) is stronger 
than the effect that eccentric mergers become more unlikely with 
increasing a\ if ei were fixed. 

As Table [3] shows, for a small fraction of MS mergers (i.e. 
about 0.009 of all MS mergers for TSM1) the merger does not oc- 
cur due to high eccentricity, i.e. the merger occurs in a circular or- 
bit. In these cases KCTF is responsible for strong orbital shrinkage 
leading to a tight inner orbit, eventually resulting in a merger in a 
circular orbit due to GW emission within a Hubble time. The frac- 
tion of MS merger systems in which this occurs (0.009 for TSM1) 
is indeed comparable with the fraction of MS merger systems that 
satisfy our KCTF criterion on the MS (0.010 for TSM1); this is 
also the case for TSM2. We find that such strong KCTF on the MS 
occurs only in low-mass inner binary systems (mi < 1.25 M Q ) in 
which t he components possess conv ective envelopes (i.e. the sce- 
narro of lFabrvcky & Tremainell2007l) . 

The eccentric MS mergers, which constitute the bulk of MS 
mergers, typically occur early on the MS, with the majority of 
mergers (about 70%) occurring within 10% of the primary MS life- 
time (Fig. [7] solid line). This is because in many cases the colli- 
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Figure 7. Cumulative distributions for TSM1 of the times t of MS merger 
(solid line) and MS destabilization (dashed line) normalized to the MS 
timescale, estimated as tms = 10(mi/M s ) Gyr. 



sion occurs during the first few Kozai cycles which typically have 
timescales of a few Myr, a small fraction of the MS lifetime of the 
primary stars considered in the triple sample. The time of merger 
is thus mainly determined by the Kozai period P K cc (P2/POP2 
and hence the orbital periods. We indeed find that the distribution 
of the times of MS mergers is very similar to the initial distribution 
of the Kozai periods. The fact that most MS mergers occur early 
on the MS makes it unlikely to observe these systems as triples. 
A similar conclusion applies to the MS destabilizations (Sect. [53} 
that happen even earlier. 

A relevant question for the MS mergers is whether the merged 
MS remnant could evolve to a blue straggler star, i.e. a star more 
luminous than stars at the MS turn-o ff point. This triple b lue str ag- 
gler scenario has been proposed by IPerets & Fabrvckvl (120091) . If 
the MS stars merge very early in the evolution then the age of the 
merged MS star would appear to be consistent with its mass and 
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TSMl 



/kctf 

MS PHG/PRGB PAGB SHG/SRGB SAGB fcs fee 

TSM2 TSMl TSM2 TSMl TSM2 TSMl TSM2 TSMl TSM2 TSMl TSM2 TSMl TSM2 TSMl TSM2 

7.743 0.018 0.012 0.242 0.233 0.147 0.212 0.029 0.025 0.001 0.001 0.408 0.463 0.607 0.561 

3.986 0.010 0.007 — — — — — — — — 0.001 0.001 0.991 0.994 

3.717 0.007 0.005 — — — — — — — — 0.001 0.001 0.994 0.996 

0.128 0.000 0.000 — — — — — — — — 0.000 0.000 1.000 1.000 

0.124 0.000 0.000 — — — — — — — — 0.000 0.000 1.000 1.000 

1.467 0.035 0.022 0.650 0.569 0.209 0.283 0.110 0.089 0.000 0.000 1.000 1.000 0.130 0.146 

0.093 0.131 0.054 0.135 0.095 — — — — — — 1.000 1.000 0.798 0.904 

0.693 0.048 0.038 0.901 0.896 — — — — — — 1.000 1.000 0.045 0.046 

0.094 0.000 0.000 0.025 0.023 — — — — — — 1.000 1.000 1.000 1.000 

0.266 0.000 0.000 0.213 0.172 0.905 0.925 — — — — 1.000 1.000 0.002 0.002 

0.146 0.000 0.000 0.260 0.240 0.795 0.818 0.289 0.294 — — 1.000 1.000 0.000 0.000 

2.290 0.024 0.014 0.496 0.423 0.420 0.534 0.039 0.028 0.003 0.003 0.893 0.924 0.098 0.072 

0.376 0.007 0.005 0.994 0.986 0.000 0.000 — — — — 1.000 1.000 0.000 0.000 

0.139 0.079 0.059 0.903 0.916 0.000 0.000 0.000 0.000 — — 1.000 1.000 0.000 0.000 

0.110 0.000 0.000 0.986 0.979 0.000 0.000 0.087 0.060 0.000 0.000 1.000 1.000 0.000 0.000 

0.282 0.106 0.067 0.188 0.276 0.198 0.291 — — — — 0.277 0.435 0.607 0.486 

0.099 0.027 0.011 0.444 0.473 0.559 0.606 0.023 0.012 — — 1.000 1.000 0.082 0.046 

0.188 0.004 0.006 0.287 0.253 0.790 0.852 0.069 0.032 — — 1.000 1.000 0.018 0.009 

0.008 0.000 0.000 0.000 0.000 0.132 0.153 0.012 0.010 — — 1.000 1.000 1.000 1.000 

0.178 0.000 0.000 0.127 0.090 0.918 0.932 0.063 0.042 0.000 0.000 1.000 1.000 0.000 0.000 

0.746 0.001 0.000 0.257 0.202 0.793 0.845 0.063 0.038 0.000 0.000 1.000 1.000 0.000 0.000 

0.128 0.000 0.000 0.130 0.104 0.719 0.786 0.145 0.105 0.025 0.022 0.873 0.899 0.121 0.091 

81.400 0.002 0.001 0.067 0.046 0.205 0.191 0.015 0.011 0.027 0.026 0.047 0.046 — — 

75.836 0.000 0.000 0.029 0.015 0.178 0.162 0.010 0.007 0.024 0.022 0.000 0.000 — — 

5.104 0.005 0.003 0.536 0.470 0.562 0.614 0.080 0.071 0.081 0.090 0.611 0.655 — — 

0.460 0.150 0.096 0.510 0.490 0.330 0.416 0.033 0.025 0.001 0.001 0.794 0.855 — — 

10.857 0.000 0.000 0.001 0.001 0.024 0.030 0.000 0.000 0.001 0.001 0.000 0.000 — — 

3.722 0.000 0.000 — — — — — — — — 0.000 0.000 — — 

5.124 0.000 0.000 0.001 0.001 0.027 0.031 0.000 0.000 0.000 0.000 0.000 0.000 — — 

4.603 0.000 0.000 0.001 0.001 0.030 0.034 — — — — 0.000 0.000 — — 

2.011 0.000 0.000 0.002 0.001 0.064 0.085 0.002 0.002 0.003 0.004 0.000 0.000 — — 

0.769 0.000 0.000 0.002 0.002 0.036 0.049 — — — — 0.000 0.000 — — 

0.977 0.000 0.000 0.002 0.000 0.097 0.115 0.004 0.002 0.007 0.007 0.000 0.000 — — 

0.122 0.000 0.000 0.002 0.001 0.066 0.092 0.001 0.002 0.009 0.005 0.000 0.000 — — 



Merger 8.054 
• On MS -> single MS 4.547 



CO WD 

■ ONe WD 

■ CC SN -> NS 



4.377 
0.071 
0.064 



■ Post-MS 



1.379 



• HG + MS (2) 0.074 
. RGB + MS (3) 0.755 

• CHeB + MS (4) 0.082 

• AGB + MS(5) 0.177 

• AGB + CHeB (4) 0.112 

> Compact object 2.128 



i He WD 
• He WD 
■ He WD 
i CO WD 
i CO WD 
>COWD 
i CO WD 
>COWD 
i CO WD 
' CO WD 



+ MS (3) 
+ RGB (7) 
+ He WD (7) 
+ MS (5) 
+ HG (5) 
+ RGB (5) 
+ CHeB (4) 
+ HeMS (8) 
+ He WD (8) 
+ CO WD 



0.471 
0.156 
0.126 
0.297 
0.060 
0.149 
0.007 
0.124 
0.618 
0.095 



No merger 



81.717 



> fli,//AU > 12 75.535 

> 10~ 2 < ai,//AU < 12 5.614 

> ai,//AU < 10~ 2 0.569 



Triple destabilization 



10.228 



. On MS 
> Post-MS 

• AGB + MS 



3.680 
4.635 

4.036 



> Compact Object 



1.913 



i CO WD + MS 
. CO WD + AGB 
> CO + CO WD 



0.888 
0.751 
0.122 



Table 3. Probabilities P (in per cent) of various evolutionary channels for both sampled triple populations. Any two objects in the first column refer to inner 
binary components. For the post-MS and compact object mergers the integer numbers in brackets denote the stellar type of the remnant according to the 
binary evolution algorithm that occurs most frequently. We refer to Table |2]for the meaning of the stellar types and for abbreviations used and to Fig. |5]for 
a complete overview of all occurring combinations of stellar types prior to an inner binary merger. In the no merger channels a\j denotes a\ at t = In- The 
quantities /kcTF: fcE and /ec denote the fraction of systems within each channel for which significant KCTF, CE evolution and a merger as a result of a highly 
eccentric collision applies respectively. For /kctf a distinction is made between the type of the star in which significant tidal energy is dissipated: MS (both 
stars are MS stars or the primary is a compact object and the secondary is a (low-mass) MS star), PHG/PRGB/PAGB (primary star is a HG/RGB/AGB star) 
and SHG/SRGB/SAGB (secondary star is a HG/RGB/AGB star). A long dash ( — ) indicates that this particular combination is not applicable. 



so it would not be observed as a blue straggler star. For this rea- 
son most of the MS mergers in our sample are unlikely to be blue 
straggler progenitors. 

5. 1.2 Post-MS mergers 

For any merger after the MS an evolutionary change in the inner 
binary system plays a key role in the mechanism that leads to the 
merger. We progressively discuss the various post-MS merger sce- 
narios of Table [3] The radius expansion during the primary HG is 
responsible for an eccentric collision for most mergers during this 
phase. Note that in such an event the Binary _c algorithm invokes 



CE evolution, in contrast to the eccentric MS mergers. For about 
20% of the HG + MS mergers strong KCTF triggered by the sub- 
stantial increase of the tidal strength quantity k am j/T t during the 
primary HG phase leads to orbital circularization and shrinkage to 
«[ « 0.4 AU. Subsequently, CE evolution is invoked, leaving a sin- 
gle HG star. For the RGB + MS mergers one might expect eccen- 
tric collisions to be important because the RGB phase is associated 
with a substantial increase in radius. However, the large increase in 
the tidal strength quantity k dm j/Ti during this phase is responsible 
for strong KCTF which circularizes and shrinks the inner orbit to 
a t « 0.3 AU leading to CE evolution, which happens in nearly all 
cases. 
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During the CHeB phase (radiative envelope) tides are unim- 
portant and most mergers are driven by radius expansion in con- 
junction with eccentricity cycles, in all cases leading to an eccen- 
tric merger. As expected, KCTF is once again important for pri- 
mary AGB mergers (convective envelope). For these systems the 



secondary at the moment of merger is most likely a MS or a CHeB 
star. For primary AGB mergers KCTF is effective mainly during the 
primary AGB phase, although in about 25% of these systems KCTF 
is also important during the primary RGB phase. After KCTF a CE 
is invoked; because the orbit at this point is tight the inner binary 
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Figure 8. Distribution for TSM1 of /J = 02/^1 prior (solid line) and after 
(dot-dashed line) the first CE for all systems in which CE evolution occurs 
and that survive the CE phase. Shown with vertical lines are the median 
values (solid) and spread values (dashed) of /J cl it,GR (Eq.[4j for the systems 
in which CE occurs, just after the CE (96% of all cases fall within the in- 
dicated spread). Lastly, the dotted line shows the distribution of fi prior to 
merger for all eccentric compact object mergers. 



does not survive the CE phase and merges into an AGB star (in 
case of AGB + MS merger) or CHeB star (in case of AGB + CHeB 
merger). 

The latter AGB + CHeB mergers, as well as RGB + HG and 
RGB + RGB mergers, have initial distributions of q x peaked to- 
wards high values (q t > 0.9), explaining the sharp increase of 
/post-Ms merge for 0.9 < q\ < 1 (Fig. In addition, Fig. [6] shows 
that post-MS mergers typically have small /3 similarly to MS merg- 
ers, although /post-Ms merge shows a clear tail for larger ft which is 
not the case for /ms merge; for compact object mergers this tail is 
even more prominent. Furthermore, post-MS mergers tend to have 
smaller initial a\ than the MS mergers. This can be understood by 
noting that KCTF is typically more effective for tighter inner binary 
systems. 



5.1.3 Compact object mergers 

If KCTF triggers a CE during the inner binary giant phases and 
the inner binary system survives this CE, a possible outcome is 
a compact object merger with the primary either a He WD (KCTF 
during RGB) or a CO WD (KCTF during RGB/AGB). Fig.[8]shows 
the distribution of p just prior to and after the first CE for all systems 
in which CE evolution occurs and that survive the CE (this includes 
systems in which an inner binary merger does not occur). After the 
first CE p ~ 10 5 , which is much larger than the critical value /? cr j t ,GR 
(cf. Eq. f4]l for which general relativistic precession dominates in 
these systems (Fig. [8] vertical lines). We therefore find that after 
the first CE p is large enough that subsequent Kozai cycles do not 
occur. In other words, the inner binary evolution after the first CE 
is not affected by the tertiary. 

This subsequent evolution is driven by stellar evolution of the 
secondary (possibly invoking an additional CE) and/or GW emis- 
sion if the inner orbit is sufficiently tight. The type of merger is 
mainly determined by the strength of KCTF: the stronger KCTF, 
the tighter the inner orbit after the first CE and the smaller the 
probability for the secondary to become a highly evolved star. 
This results in a multitude of compact object mergers, among 
which are He WD + HG/RGB mergers and CO WD + He WD 
mergers. The former are believed to be important channels lead- 
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Figure 9. Distribution for TSM1 of the initial a\ for compact object merg- 
ers in which strong KCTF triggers a CE and finally a merger in a circular 
orbit ("circ"; solid line) and in which an eccentric collision occurs ("ecc"; 
dashed line). In addition, for the former systems the distribution is shown 
of a\ after KCTF during the primary AGB phase (dot-dashed line); for the 
latter systems the distribution is shown of a\ just prior to the eccentric colli- 
sion (dotted line). Note that the systems having a\ ~ 10 _1 AU after primary 
AGB KCTF are those in which strong KCTF occurred also during an earlier 
stage, most notably the primary RGB phase. 



ing to carbon-ric h K-type giants, known as early-type R stars 
jlzzard et al1l2007l) ; the latter have been discussed as possible pro- 
genito rs of H-deficient ca rbon-rich supergiants known as RCrB 
stars dTisserand et~aL1l2009ri . 

For a fraction f EC » 0.1 of all compact object mergers strong 
KCTF and CE evolution during the MS and post-MS phases are 
avoided. These systems are typically initially wide (ai > 10 2 AU; 
cf . Fig. [9] dashed line) and expand to even wider orbits during the 
AGB phases due to wind mass loss (Fig. [9] dotted line), as op- 
posed to shrinking as a consequence of KCTF. Similar to the ex- 
ample system in Sect. 13.21 the ratio p = ail<i\ decreases because 
of this mass loss. Fig. [8] (dotted line) shows that the resulting p 
is narrowly distributed around p as 15. Octupole order effects are 
very important in such systems, with |e oct | peaking around « 10" 2 . 
In addition, because at this stage the inner orbit semi-major axis 
is large (typically a t ss 10 3 AU; cf. Fig. [9] dotted line) there is 
no significant damping of Kozai cycles due to additional sources 
of inner orbit precession. Consequently, in such systems orbital 
flips can occur as in Sect. 13.21 durin g which extremely high ec- 
centricities can be reached (e.g.lNaoz et al .l201ll : lLithwick & Naozl 



l201ll ; IShappee & Thompsonll2012h . in our calculations as high as 
1 -ei ~ 10" 8 . 

In most cases these high eccentricities lead to orbital colli- 
sions. The value of the eccentricity required for such a collision 
is determined by «i and the radius and hence the type/mass of the 
secondary object (mainly a MS star, a CHeB star or a CO WD) and 
varies between 1 - e\ ~ 10~ 3 (CHeB companion), 1 - e\ ~ 10~ 5 
(MS companion) and \ — e\~ 10~ 7 (CO WD companion). We will 
discuss the CO WD + CO WD eccentric merger scenario in more 
detail in Sect. [6] 

In a few cases, mainly for CO WD + low-mass MS systems, 
the high eccentricities do not lead to eccentric collisions but induce 
very strong tidal friction in the low-mass secondary that possesses 
a convective envelope, thus circularizing and shrinking the inner 
orbit to a\ ~ W 2 AU. These cases correspond to the MS /kctf 
fractions for the CO WD + MS merger channel and the no merger 
channel with aij < 12 AU in Table [3] Depending on the precise 
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log 10 (ai/AU) 

Figure 10. Distribution for TSM1 of the initial (solid line) and final a\ 
(dashed line) for all no merger systems. 

values of a\ and the masses the subsequent tight binary does or 
does not merge within a Hubble time due to GW emission. The 
systems in which the latter scenario applies constitute about 20% 
of the tightest no inner binary merger systems (aij < 10~ 2 AU 
in Table [3} . They may be recognized in Table [3] as the no merger 
systems with a^j < 1(T 2 AU in which a CE does not occur. 

5.2 No inner binary mergers 

If the initial inclination angle is not close to 90° and/or the ini- 
tial p is large then the maximum eccentricity induced by Kozai cy- 
cles during the evolution may not be high enough to drive an inner 
binary merger. This is demonstrated by Fig. [6] which shows that 
/no merge increases strongly with increasing /? (/ n0 mcrgc — > 1 for suffi- 
ciently large [}) and that systems with no or moderate inclination are 
likely not to merge. Furthermore, / nom e rg e increases with decreasing 
ci[ which might be somewhat counterintuitive. This can be under- 
stood by considering that as a t decreases /? is also generally larger, 
thus decreasing the eccentricity maxima. 

Although by definition a merger does not occur for these sys- 
tems, it is possible that the inner orbit is shrunk significantly due 
to KCTF as is illustrated by Table [3] in about 6% of all systems 
the final inner orbit semi-major axis is smaller than 12 AU, which 
is the minimum possible value in the absence of the tertiary (cf. 
Sect. 14. D . In the majority of these systems the main points in 
the evolution at which KCTF is important are the primary RGB 
and AGB phases occurring in roughly equal proportion. Fig. [TOl 
compares the initial and final distributions of a\. The final distri- 
bution shows a broad peak mainly due to strong KCTF around 
a Lf ~ lfr 15 AU ~ 0.03 AU (Pi ~ 2d). For the tightest of these 
systems (a i f < 10~ 2 AU) it is also possible that very high eccen- 
tricity cycles lead to a tidal capture if the systems consists of a CO 
WD + low-mass MS star, which is the same scenario discussed in 
the last paragraph of Sect. 15.1.31 

The final distribution of a\ shows a gap around a\j ~ 1 AU. 
This gap separates systems in which KCTF is followed by a CE 
(a t f < 1 AU) and in which a CE does not occur after KCTF (aij > 
1 AU). In the latter case KCTF is responsible for an enhancement 
of systems near final values of «i ~ 10 AU. If a CE occurs then 
this affects the final configuration of the inner binary system. This 
is shown in Table [4] where the likelihoods of the most important 
inner binary configurations at the end of the evolution are shown 
for the no merger systems. A distinction is made between the three 
groups with aij > 12 AU ("wide"), 10~ 2 < aij/AV < 12 ("tight") 
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Figure 11. Distribution for TSM1 of the initial inclination angle for the no 
merger systems (solid line) and the final inclination angle for the no merger 
systems with final a\ < 12 AU (dashed line) and a\ > 12 AU (dot-dashed 
line). 



Wide Tight Very tight 

TSM1 TSM2 TSM1 TSM2 TSM1 TSM2 

CO + CO WD 0.370 0.428 0.189 0.209 0.010 0.016 

CO WD + MS 0.599 0.535 0.503 0.494 0.391 0.392 

CO + He WD 0.000 0.000 0.077 0.093 0.255 0.293 

He + He WD 0.000 0.000 0.014 0.013 0.058 0.049 

HeWD + MS 0.000 0.000 0.188 0.159 0.284 0.246 



Table 4. Likelihood of the main inner binary configurations at the end of 
the evolution for the no merger systems. These systems are divided into 
three groups, based on the final value of ay. wide (a\t > 12AU), tight 
(10~ 2 < ai,//AU < 12) and very tight (aij < 10~ 2 AU). The fractions are 
relative to the total number of systems in each group. 



and a t f < 10~ 2 AU ("very tight"). A CE does not occur in the 
wide systems and therefore these end mainly as CO WD + MS and 
CO WD + CO WD systems. In the case of a CE there is also a 
possibility for forming a He WD. The tighter systems can therefore 
also end as CO WD + He WD, He WD + He WD and He WD + 
MS systems. 

The no merger systems are interesting observationally because 
the hierarchical structure remains intact during the entire evolution 
and the process of KCTF distinctly affects the orbital inclination an- 
gle. Typically during episodes of KCTF ej gradually approaches its 
maximum value in the Kozai cycle, followed by rapid circulariza- 
tion. The inclination angle during this process tends to remain fixed 
at its value associated with the eccentricity maximum, i.e. close to 
one of the two critical values determined by cos 2 (/ tot ) = 3/5. There- 
fore an observational marker of systems in which KCTF has once 
been im portant is an inclination angle close to these critical values 
(see also lFabrvckv & Tremainel2007|) . Fig.QT]shows the initial and 
final inclination angles for the no merger systems, distinguishing 
between a^j > 12 AU and a\j < 12 AU. The initial distribution 
shows a clear lack of systems near i t01 as 90° as also illustrated by 
Fig. [6] Systems for which strong KCTF applies at some point in 
the evolution (a lt f < 12 AU) end with inclination angles strongly 
peaked towards the critical values, whereas for systems in which 
KCTF is weaker or does not act at all (final a\ > 12AU) these 
peaks are much less prominent. We therefore expect triples with 
tight inner orbits to show a different inclination distribution than 
those with much wider inner orbits. 
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5.3 Triple destabilizations 

If the initial /3 is particularly low, typically ft < 10 (cf. Fig. |6j, 
then the triple system may at some point in the evolution become 
dynamically unstable. Subsequent evolution is not modeled by our 
triple evolution algorithm because the secular evolution equations 
no longer apply in this case and three-body simulations are needed 
for an accurate description of the d ynamical evolution. Th is sce- 
nario has been studied in detail by |Perets~& Kratter (2012), who 
refer to it as the triple evolution dynamical instability. In our sim- 
ulations the destabilization occurs mainly during the MS or when 
the primary is an AGB star or a CO WD (cf. Table. [3}. 

In the case of MS destabilization the triple system is initially 
marginally stable (i.e. only just satisfies yS > p a \ t , cf. Sect. |2.3l l but 
due to octupole-order terms of the STD, which are important since 
P is (very) small and/or e 2 is high. e 2 varies periodically until it 
reaches a value high enough such that [} < /? crit , i.e. a triple desta- 
bilization. The time when this occurs is determined by the Kozai 
period P K . Similarly to the MS mergers, this occurs early in the 
evolution with most (90%) destabilizations occurring within 10% 
of the primary MS lifetime (cf. Fig. 0. 

In the other cases destabilization is triggered by mass loss in 
the inner orbit which, if fast and isotropic, acts to decrease /? (i.e. 
the same mechanism discussed in the context of eccentric compact 
object mergers in Sect. l5.lt to a point where fj < This happens 
when the primary loses a significant amount of mass as it evolves 
from the AGB phase to a WD and similarly when this happens to 
the secondary. In a small number of cases both inner binary compo- 
nents are CO WDs when the instability occurs and since there ex- 
ists a finite probability of collision in the tri ple evolution dynamica l 
instability (approximately 0. 1 as found bv lPerets & Kratterll2012h 
this could potentially lead to a CO WD collision. This scenario is 
included in Sect. [6] 
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Figure 12. Distribution for TSM1 of the inner binary total mass mm = 
mi +»i2 (top; vertical dashed line denotes Mq,) and mass ratio qi = m2ltn\ 
(bottom) just prior to CO WD merger or triple destabilization, for the CO 
WD inner binary circular (solid line) and eccentric (dashed line) mergers 
and the CO WD destabilizations (dotted line). 



6 IMPLICATIONS OF TRIPLE-INDUCED CO WD 
MERGERS 

As mentioned in the introduction the CO WD merger channel is 
con sidered as an important can did ate progenitor chann el for SNe la 
(see lMaoz & Mannucci|[20ll and lWang & Har]|2012l for recent re- 
views). In this section we explore the implications of triple-induced 
CO WD mergers by estimating the expected SN la rate. First we 
briefly discuss the triple-induced CO WD merger scenarios that we 
find in our population synthesis study (Sect. 16. W . In order to esti- 
mate the SNe la rate we must make assumptions for which com- 
binations of CO WD total masses and mass ratios prior to merger 
lead to a SN la explosion. We base these assumptions on the re- 
sults of various detailed simulations of CO WD mergers (Sect. [6T2l . 
With these assumptions and by normalizing the results we obtain 
the expected SNe la rates. We then compare these rates to the bi- 
nary population synthesis study of Claeys et al. (in prep.) and to 
those inferred from observations (Sect, |6~3~| . 

6.1 Scenarios 

In Sect. [5] we have found that in the sampled triple populations 
CO WD mergers can be the result of 1) strong KCTF followed by 
one or multiple CE phases eventually followed by a merger in a 
circular orbit due to GW emission, 2) extremely high eccentric- 
ities (1 - ei ~ 10~ 7 ) excited in wide inner binaries with a distant 
companion where KCTF is avoided, resulting in an orbital collision 



and 3) triple destabilization driven by mass loss du ring the inner bi- 
nary A GB phases, possibly leading to a collision jPerets & Kratter] 
2012). Note that systems in the first scenario typically have initial 
a,\ in the range 10 1 < ai/AU < 10 2 whereas for the second sce- 
nario this range is 10 2 < ai/AU < 10 3 (cf. Fig.|9j. As shown in 
Table [3] scenario 1 is about 10 times more likely than scenario 2 
(i.e. / EC k 0.1 for the CO WD merger channel), whereas scenario 
3, taking into account that a triple instability leads to a collision in 
about 0. 1 of all cases jPerets & Krattej2012l) . is about 10 times less 
likely than scenarios 1 and 2 combined. Note that in the absence of 
the tertiary these channels would not exist for the sampled popula- 
tions, i.e. they are all induced by the secular gravitational influence 
of the tertiary. In binary population synthesis studies such channels 
are therefore not taken into account. 

6.2 Likelihood of a SN la explosion 

Many uncertainties exist as to which configurations of CO WD 
mergers lead to a SN la explosion. Nevertheless, estimates of the 
minimum inner binary total mass m m = m\ + m 2 and mass ratio 
q { = m 2 /mi required can be inferred from detailed simulations. 
Here it is important to distinguish between mergers occurring in 
circular orbits (with GW emission driving the merger) and mergers 
occurring in highly eccentric orbits, very similar to head-on col- 
lisions. The latter are more likely to result in a SN la explosion 
because of the high relative speed, on the order of the escape speed 
(~ 10 3 km/s), which causes shocks that aid the ignition of the ther- 
monuclear explosion. For circular mergers the minimum q\ could 
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be approximately q\ as 0.8 (e.g. lPakmor et al ] l201ll) provided that 
m m > M a ~ 1.4 M . Furthermore, for eccentric mergers a min- 
imum mass ra t io of q\ = 2/3 can be inferred fro m Table 1 of 
iRosswog et al.l d2009h and qj = 0.6 from Table 2 of Rask in et al.l 
1 201 Oh . According to the latter study the minimum total mass is 
between 1.0 and 1.2 M Q . 

Fig [12] shows the distributions of m m (top) and q t (bottom) 
prior to merger for the three triple-induced CO WD merger sce- 
narios. In the circular merger scenario the merger can be both sub- 
Mo, (57%) and super-M C i, (43%) and q^ is broadly distributed with 
0.4 < q\ < 1.0. In the eccentric merger scenario mm exceeds the 
upper limit of the minimum m ass for a SN la event of 1.2 M Q ac- 
cording to Rask in et al. for most systems (95%), whereas qi 
is peaked towards high values, 0.6 < qi < 1. Lastly, in the destabi- 
lization scenario m 1B is typically slighter lower than in the eccentric 
merger scenario whereas qi is even more strongly peaked towards 
high values. 

The fact that qi in scenario 1 can be moderately small suggests 
that not all mergers in this scenario lead to a SN la explosion. In 
binary population synthesis st udies a distinction in mass ratio is 
not usually made, however (see lChen et al ] |2012l for a study of the 
effect on the binary population synthesis rates of including more 
detailed constraints). In order to facilitate a comparison with results 
from such binary population synthesis studies we shall therefore 
disregard the mass ratio criterion and only consider a restriction on 
mm '■ in Sect. l6.3l we shall assume that a SN la explosion is the result 
in scenario 1 (circular merger) if m m > M cb and in scenarios 2 and 
3 (collision) if mm > 1.2 M e . 

6.3 Expected SNe la rates 

Conforming to the procedure commonly employed in binary pop- 
ulation synthesis studies, we normalize the number of systems in 
which a SN la event is expected to occur to the part of the galactic 
mass represented by the sampled systems. The details of this nor- 
malization are included in Appendix lAl here we mention our most 
important assumptions. Firstly, it must be taken into account that 
the sampled triple populations constitute only a small fraction of all 
triple systems (cf. Sect. 14.1) . In addition we must assume specific 
values for the binary and triple fractions of the galactic population. 
The binary fraction is a function of mass and is found to ran ge 
between ff bin = 0.56 for s olar-like stars jRaghavan et al.ll2010l) to 
a bin > 0.7 for O & B stars dKouwenhoven et al J 20071) . For simplic- 
ity we assume here an intermediate constant value of = 0.60. 
Statistics of the triple fraction a lT , which is similarly a strong func- 
tion of mass, are even less well-constraine d. The observed triple 
fractions range fr om 0. 1 1 for solar-like stars (Raghava n et alj20ld : 
iRastegaevI 2010) to 0.5 for more massive B stars (Evans 201 II) . 
Again neglecting the mass dependence we take a conservative inter- 
mediate constant value a tl = 0.25, where a tl is not to be understood 
as a subset of a bin . We neglect higher-order multiplicities, thus giv- 
ing a single star fraction of 1 -a bin -a tr = 0. 15. With these fractions 
we calculate the total mass represented by the number of triple sys- 
tems we have calculated and with which we normalize the SNe la 
rates for TSM1 (Eq. |A5) . We express the result in rate per 10 10 M Q 
per century (SNuM). Fig.ll3lshows these rates as a function of time 
t after an assumed starburst at t = 0, i.e. the delay time distribution 
(DTD). The contributions of the three scenarios of Sect. |6.1| taking 
into account the assumptions of Sect. l6.2l are indicated separately. 

Fig.[T3]shows that the main contribution to the triple-induced 
rate is from the inner binary circular merger channel. The inner bi- 
nary eccentric merger channel contributes significantly at late times 
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Figure 13. SNe la delay time distribution (DTD) for TSM1 of the triple- 
induced inner binary CO WD merger and destabilization channels with the 
criteria for producing a SNe la event discussed in Sect. 16.21 Shown are 
the combined rate (heavy solid line), the rate due to inner binary circu- 
lar mergers with mm > Mai (thin solid line) and eccentric mergers with 
mm > 1 -2 M B (dashed line) and the rate due to CO WD collisions induced 
by triple destabilizations assuming that in 10% of all inner binary CO WD 
destabilizations a collision subsequently occurs (dotted line). 



(r ~ lOGyr), whereas it does not at earlier times. This is because 
after the CO WD formation the Kozai period P K is typically quite 
long (the distribution of P K peaks around P K ss 0.3 Gyr; see also 
Sect. 13. 21 and it typically takes at least a few cycles before an or- 
bital flip occurs (cf. Fig. El bottom right panel). Such an orbital flip 
plays a crucial role in driving eccentricities high enough for a col- 
lision in these systems (1 - e\ < 10~ 7 ). The destabilization channel 
does not contribute significantly to the total rate at any time. 

Furthermore, we show in Fig. [14] the combined expected 
triple-induced SNe la rate according to the three channels discussed 
above (dashed line), the predicted SNe la DTD from the binary 
population synthesis st udy of Claeys et a l. (in prep.) (solid line) 
and the observed rates dMaoz et al.[|2012l) . Note that in Claeys et 
al. (in prep.) a binary fraction of a b [ n = 1 is assumed and that con- 
tributions are shown from both the single degenerate channel (sin- 
gle CO WD accreting material from a non-degenerate donor star 
until it reaches Mqu and explodes; not taken into account for the 
triple population here) and the double degenerate channel (circu- 
lar CO WD binary merger due to GW emission with % > M C h)- 
Fig. [13 shows that the triple channel DTD generally follows the 
same shape as the binary population synthesis DTD, with a rate 
proportional to t for times later than t x 0.4 Gyr (the best fit- 
ting slope for the triple case for t > 0.4 Gyr is w - 1 .08). In absolute 
terms, however, the triple-induced rates are much lower than the bi- 
nary population synthesis rates: the former are typically factors of 
~ 10 2 - 10 3 lower than the binary population synthesis rates. The 
total time-integrated rates are about 4 • 10~ 4 M G ' and 2 • 10~ 6 M Q ' 
for the binary and triple cases respec tively while the observed in- 
tegrated rate from lMaoz et al.l d2012l) is (1.30 ± 0.15) • 10~ 3 
The triple-induced channels discussed here are therefore not able 
to resolve the dis crepancy between the pre dict ed and observed SN e 
la rates (see also lMaoz & Mannuccill2012l and lWang & Hanll2012l) . 

It is important to emphasize, however, that the present study 
does not give a complete picture of the triple-induced rate because 
we have not taken into account triple systems with inner binary or- 
bits with li < 12 AU. Further study is therefore needed to examine 
the contributions from triple systems with l\ < 12 AU. 
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Figure 14. SNe la DTD according to all triple-induced channels of Fig. 
1 131 (dotted line) and according to the binary population synthesis study of 
Claeys et al. (in prep.) (solid line). In addition, observed SNe la rates (points 
with error bars) are shown from Maoz et al. (2012). 



7 DISCUSSION 

7.1 Accelerating CO WD mergers through Kozai cycles 

In Sect.[T]we mentioned the scenario o f lThompsonl l l20T ]]) in which 
the tertiary induces high-amplitude Kozai cycles in close CO WD 
binary systems, thus reducing the GW merger time and producing 
a GW signature that is distinct from th e signature from GW merg- 
ers in circular orbits. Thom pson! d201 lh considers triples with close 
inner binary CO WD systems (10~ 2 < ai/AU < 10) and the outer 
orbit is assumed to be fairly tight (10 < a%ja\ < 100). For coeval 
triple systems it is also important to take into account the evolution 
prior to the formation of the binary CO WD system. Such evolu- 
tion likely involves CE evolution in the inner binary system, which 
significantly shrinks the inner orbit and widens the outer orbit, thus 
substantially increasing the semi-major axis ratio [} = ai/ai. For 
the coeval triples considered in this work we find that the semi- 
major axis ratio /? is large enough (J3 ~ 10 5 , cf. Fig.[8]( to completely 
damp any subsequent Kozai cycles because of general relativistic 
orbit precession in all systems in which a close CO WD system 
is form ed after a CE. We thus find that the scenario of Thompson 
d201ll) does not occur in coeval triples with l { > 12 AU (the sys- 
tems with 1 1 < 12 AU have not been modeled here). 

Furthermore w e briefly mention the very recent work of 
iKatz & Pond 420121") who find that collisions of compact objects 
can occur in triple systems in which high eccentricities are reached 
through Kozai cycles. This scenario is very similar t o the eccentric 
merge r scenario found in this study (cf . Sect. 15.1.3} . Katz & Dong 
( 120121) use direct iV-body methods to model the dynamical evolu- 
tion of triple WD systems and find that the secular description of 
Kozai cy cles breaks down in m any of their systems. We note, how- 
ever, that lKatz & Donal l2012h only consider equal-mass inner bi- 
naries, whereas we find that most eccentric WD collisions occur in 
systems with qi < 1 (Fig. [T2] bottom panel) in whi ch octupole- 
order terms dominate the dynamics. In addition, IKatz & Pond 
<2012h do not take into account stellar, binary and triple evolution 
prior to the formation of the double WP system. We find that these 
aspects are very important, however: most systems with high a\, 
low /3 = a%ja\ and high inclination, i.e. potential progenitors for 
the collision scenario, experience an early merger already during 
the MS prior to evolving to systems with compact objects (cf. Sect. 

EED. 



7.2 Uncertainties in the triple evolution algorithm 

In the triple algorithm developed in this work (Sect. [2) various sim- 
plifying assumptions are made, some of them not generally valid. 
Most importantly, we assumed a fast and isotropic wind mass loss 
for the outer binary orbit. If the inner orbit is tight and the inner 
orbital speeds are fast compared to the stellar wind speeds then the 
wind from either component in the inner binar y system is likely 
affect ed by its inner binary companion (e.g. Ide Val-Borro et all 
2009), thus invalidating our assumption. In our triple population 
this could be the case for some of the tighter inner binary systems 
containing AGB stars (with typically slow winds), in which case 
strong wind interaction is possible even in orbits with semi- major 
axes exceeding 12 AU (e.g. Mohamed & Podsiadlowski 2007). Be- 
cause our results are strongly dependent on this assumption of fast 
and isotropic wind mass loss it would be worthwhile to implement 
a more sophisticated treatment of mass loss in our triple algorithm. 
A similar problem arises whenever a CE is invoked: in the present 
study it is assumed that the CE material from the inner binary is lost 
as if it were a fast and isotropic wind from the inner binary system. 

As Fig. [6] shows our results are strongly dependent on the as- 
sumed triple distributions which are generally poorly constrained. 
Although the results from the two populations TSM1 and TSM2 are 
very similar (Table|3}, it must be noted that in both populations the 
initial inclination angle is sampled from the same distribution (i.e. 
a uniform distribution in cos[i tol ]). More observational constraints 
are needed on the initial triple parameters, in particular f tot . 

Furthermore, we assumed that the secular three-body dynam- 
ics, tidal friction and GW emission can be described in concert 
by linearly adding the relevant OPE terms (Eq. |3}- It would be 
of value to investigate the validity of this assumption by means of 
more detailed simulations where a direct three-body code is cou- 
pled with hydrodynamic s and/or stellar evolution co des, e.g. with 
the AMUSE framework jPortegies Zwart et al.ll2012l) . 

Another point of concern is whether our secular approach to 
hierarchical triple systems is still valid in the limit of the extreme 
eccentricities found in the population synthesis study (for eccentric 
CO WO mergers typically \-e\ ~ 10~ 7 , cf. Sects.|3T2land |5~T3l To 
investigate this we have computed the gravitational dynamics for a 
system in our population synthesis sample after the formation of 
the triple CO WO system with a high-precision direct iV-body code 
(Boekholt, in prep.) and we indeed find such high eccentricities. 

In addition, we note that uncertainties remain in the physics 
of tidal friction, in particular in the dynamical tide model where 
high eccentricities may induce co mplicated coupling of stel l ar os- 
cillati ons with the tidal potential dWitte & Savoniiel ll999: Zahn 
2005), thu s possibly making t he tidal strength quantities fc am ,,/r, 
I Hu7ll98ll : iHurlev et alj 12002) dependent on orbital eccentricity. 
We have investigated possible consequences of this scenario by 
means of an ad hoc approach, increasing the tidal strength quan- 
tities k am j/Tj by a factor of 10 3 separately for the cases of radiative 
and degenerate damping and rerunning the calculations. We find 
that this makes no significant difference. 

Lastly, a possible concern is the fact that in our treatment of 
the secular dynamics of hierarchical triple systems not all post- 
Newtonian (PN) terms of the lowest order, i.e. the 1PN terms 
<x (v/c) 2 , are incl uded in the Hamiltonian. As shown recently by 
iNaoz e t al. ( 2012), a correct derivation starting from the three-body 
Hamiltonian accurate to 1PN order introduces an additional inter- 
action term in the secularly averaged Hamiltonian. In an additional 
run of the calculations we have taken this interaction term into ac- 
count and we find no significant difference. 
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8 CONCLUSIONS 

We have performed a population synthesis study of triples with in- 
ner binary systems having initially l t = a\ ( 1 - e\\ > 12 AU and 
a primary mass range 1.0 - 6.5 M Q . Our main conclusions are as 
follows. 

1 . In the considered triple sample the inner binary system does not 
interact in the absence of the tertiary. However, with the tertiary 
present the inner binary evolution is significantly affected in about 
24% of all simulated systems. As expected the process of Kozai 
cycles with tidal friction (KCTF) can be very important during the 
primary giant phases (RGB/AGB) whereas this was not the case 
during the main sequence. KCTF during the giant phases can drive 
strong orbital shrinkage and is often followed by common envelope 
evolution. As a consequence many different merger scenarios are 
possible. 

2. After a common envelope driven by KCTF and radius expansion 
further Kozai cycles are suppressed by general relativistic preces- 
sion in the inner orbit. For close white dwarf binary systems formed 
in this scenario, this excludes the possibility of significantly reduc- 
ing the gravitational wave merg er time owing to high-amplitude 
Kozai cycles, as suggested bv lThompsonl J201 ll ). In the scenario of 
the present work the evolution prior to the white dwarf formation 
is strongly influenced by the tertiary, however. 

3. In some systems with wide inner binaries KCTF is avoided dur- 
ing the giant phases, but a subsequent decrease in the semi-major 
axis ratio a 2 /a t due to mass loss drives Kozai cycles with extreme 
eccentricity amplitudes. A subset of these systems is expected to 
experience an orbital collision between two carbon-oxygen white 
dwarfs in the inner binary system, which is a novel candidate SNe 
la progenitor scenario not found in any binary evolution path. 

4. We have identified three triple-induced carbon-oxygen white 
dwarf merger channels in our population synthesis study, among 
which the novel eccentric collision scenario, and have estimated 
the expected SNe la rates resulting from such mergers as function 
of time. We have compared these results to the binary population 
synthesis study of Claeys et al. (in prep.) and to the observations 
and find that the contribution from the triple-induced channels is 
small. We point out that further study is necessary in which triples 
with /[ < 12 AU are also taken into account. However, in this work 
we have determined a lower limit to the contribution from triples to 
the SNe la rate. 
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APPENDIX A: DTD NORMALIZATION CALCULATIONS 

We describe our method to normalize the delay time distribution 
(DTD) of Sect. |6.3| to the part of the mass represented by the sam- 
pled triple systems. Note that for the purposes of a comparison of 
the triple-induced rates to those of the binary population synthesis 
study of Claeys et al. (in prep.) a restriction must be made to the 
TSM1 population which is based on similar (inner) binary distribu- 
tions. In this work the triple populations are split into two parts with 
1.0 < m,/M < 2.0 (TSM1A) and 2.0 < m,/M < 6.5 (TSM1B) 
(cf. Sect. l4.2t . Therefore we perform the calculations separately for 
these separate mass ranges; subsequently we add the contributions 
from both mass ranges. 

We decompose the galactic population of N tot gravitationally 
bound stellar systems into A'bin = ffbinMot binary star systems, 
N lr = a? tT N tol triple star systems and N s = (1 - a'bin ~~ ttttWot single 
star systems (i.e. we neglect higher-order multiplicities). We disre- 
gard from the calculated triple systems in TSM1 the MS mergers 
and MS destabilizations because these are not expected to be part 
of the observed population. There thus remain A^ ca i c- A = 915,395 
and A^i^b = 923, 660 computed systems for the two mass ranges 
respectively. These remaining triple systems constitute a fraction 
/ ca i c of all N u triple systems such that iVeaic = / ca i c Af tr . Assuming 
that the initial primary mass mi is uncorrected with a\ and e\, 
/oalc = /cak. m , x / ca ic. c ,. C2 , ni , n2 , where f^ mi is the fraction of cal- 
culated systems to all systems with respect to the mass distribution 
and /ciiicc a j,n 2 ls this fraction with respect to the eccentricity and 
semi-major a xes distributions of the inner and outer orbits. For a 
iKroupa etal]ll993l) initial mass distribution that is used in TSM1, 

dN ( m" 1 , m Kl l < m < m Kl , 2 ; 

— oc \ m" 2 , m K r,2 < rn < m Kr ,3; (Al) 
am „, 

^ m 3 , OTKr,3 < m < mKr,4, 

where a-, = {-1.3, -2.2, -2.7) and m KlJ /M = (0.1,0.5, 1, 80), the 
fraction of systems with m\j <m\ < mi jU is given by: 



/caIo,m, = — 7 C Kr C K r,3 ("CT' ~ m Tl +> ) • ( A2 ) 

where 

c * {l^nK.-. -<•)} < A3 > 

and Cicr,; = {^>m^~" 2 ,m^T£ 2 m Kr3" 3 l- Here m\ ; = 1.0 M Q and 
m Uu = 2.0 M Q for TSM1A and m u = 2.0 M s and m hu = 6.5 M Q for 
TSM1B. Evaluating Eg. IA2I numerically for the two mass ranges 
A and B we find / cak , m ,, A ~ 6.34 • 10~ 2 and / ca ic.„„.B ~ 2.44 ■ 10~ 2 . 
Note that this implies relative weights of 0.722 and 0.278 for the 
two mass ranges respectively; these weights are used for the calcu- 
lations in Sect. [5] 

The determination of / C aic,ei,e,,ai,a 2 ls complicated by the fact 
th at a\ and a 2 are corre l ated because the stability criterion is used 
of iMardling & Aarsethl fcOQll) . This criterion has a non-trivial de- 
pendence on e 2 and q 2 = m^j(m\ + m 2 ). In addition, the selection 
l\ = a\ ( 1 - e 2 ) > lij = 12 AU is made. In an analytic treatment 
one would therefore have to integrate the TSM1 probability den- 
sity function associated with /c a ic, ei ,e2,ai,a 2 , d 4 N/ (deide^daid^) K 
e\eil(fl\Oq), with respect to e,, e 2 , «i and a 2 , where aiy8 cr j t < a 2 < 
a„, ail [\ - efj < a l < a„, < e\ < (1 - hj/a u ) 112 and < e% < 1. 
Here a t = 5R Q and a„ = 5 ■ 10 6 R a are the lowe r and upper limits re- 
spect ively for the semi-major axis distribution dKouwenhoven et ail 
2007); the upper limit of ey is a consequence of the requirement that 
a u (l - e 2 ) > In or equivalently e\ < e lM = (1 - h.ila u ) 111 . The re- 
quired integration is complicated to do analytically. 

Instead of this analytic method we determine the fraction 
,7cafc,ei,<!2,ai,a 2 by the following method: we generate a large sample 
of systems according to the method appropriate for TSM1A and 
TSM1B without making the selection h > From these systems 
we subsequently select the systems that satisfy /] > l t j. The num- 
ber of the latter systems divided by the number of sampled systems 
represents an estimate of /c a k, e] ,c 2 ,a 1 ,a, if a sufficiently large num- 
ber of systems is sampled. We perform this procedure with sizes 
varying between 10 4 and 10 6 and thus find / ca i c , e , , C2 , fl] , fl2 ~ 0. 130 for 
both TSM1A and TSM1B (i.e. this fraction is independent of pri- 
mary mass). This implies calculated fractions of / ca i C A ~ 8.24- 10~ 3 
and/ calc , B «3.17-10- 3 . 

The last step is to add the mass contributions from the single, 
binary and triple star systems. In general, for a Kroupa IMF the 
total mass M in a population of size N is given by: 

r»>K,A dN 3 c 
M = m — dm = iVx Cm > — \mJ il - mJ . 

= NMkt, (A4) 

where we have defined a "Kroupa mass" M Kl w 0.5006 M Q that 
represents the average mass of a single star in the population. Hence 
the total mass of the single star population is given approximately 
by M s ~ M Kl N s . For the binary population we take a flat mass ratio 
distribution (consistent with Claeys et al., in prep.) such that the 
secondary is on average half as massive as the primary. Therefore 
the total mass contained in the binary population, M hin , is given 
approximately by M bin ~ (l + jjM Kl N bm = fM Kr iV bin . Similarly, 
using that in TSM1 the distribution of the outer orbit mass ratio 
<?2 = m?, I {mi + mi) is flat the total mass contained in the triple 
population is M„. as (f + | ■ §)M Kl JV tl = |M Kl 7V tr . Thus the total 
mass M tot = M s + M bm + M tr represented by the number of triple 
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systems we have calculated (N all . = f ca i c N tI ) is given by: 



With a binary fraction of a bin = 0.60 and a triple fraction of 
a tr = 0.25 we thus find that the masses represented by TSM1A 
and TSM1B are given by M touA w 3.59 ■ 10 8 M and M tol>B a 
9.41 • 10 s M . We use these masses to normalize the triple-induced 
DTDs in Sect.ROl 




(A5) 



